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The spectrum of tropospheric delay fluctuations expected for a DSN antenna at 
time scales <100 s has been calculated. A new feature included in these calculations 
is the effect of aperture averaging, which causes a reduction in delay fluctuations 
on time scales less than the antenna wind speed crossing time, &5-10 s. On time 
scales less than a few seconds, the Allan deviation <j y (At) oc (A t) +1 , rather than 
<jy(At) oc (At) -1 / 6 without aperture averaging. Due to thermal radiometer noise, 
calibration of tropospheric delay fluctuations with water vapor radiometers will not 
be possible on time scales less than ~10 s. However, the tropospheric fluctuation 
level will be small enough that radio science measurements with a spacecraft on 
time scales less than a few seconds will be limited by the stability of frequency 
standards and/or other nontropospheric effects. 


L Introduction 

Radio science experiments with a spacecraft use a one- or two-way radio link between a DSN antenna 
and the spacecraft. The received amplitude and phase are used to measure quantities such as gravitational 
radiation, refraetivity profiles of the atmosphere of the planet or planetary satellite near the spacecraft, 
or the structure of a planetary ring system. Any perturbations on the link that are caused by the media 
between the DSN antenna and the spacecraft will corrupt the accuracy of the radio science measurement 
(unless the purpose of the experiment is to study that medium, e.g., the solar plasma). As the link 
frequency used for radio science experiments has increased, the relative magnitude of tropospheric and 
solar plasma phase fluctuations has changed dramatically, because plasma is dispersive at microwave 
frequencies and the troposphere is not. At S-band (2.3 GHz), solar plasma phase fluctuations dominate, 
while at Ka-band (32 GHz), tropospheric phase fluctuations will dominate except at small sun-spacecraft 
angular separations. 

Knowledge of the tropospheric fluctuation spectrum will be useful in planning and analyzing radio 
science experiments. Comparison of fluctuation levels measured by. very long baseline interferometry 
(VLBI) and water vapor radiometers (WVRs) has demonstrated that tropospheric delay/phase fluctua- 
tions at microwave frequencies are dominated by fluctuations in water vapor density [l]. 1 The spectrum 


1 C. D. Edwards, “Water Vapor Radiometer Line-of-Sight Calibration Capabilities,” JPL Interoffice Memorandum 
335.1-90-015 (internal document), Jet Propulsion Laboratory, Pasadena, California, March 30, 1990. 
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of tropospheric fluctuations has been measured by VLBI on time scales >20 s [2] and by WVRs on time 
scales >200 s [3]. A model for these fluctuations has been developed [4]. Previous calculations using this 
model (e.g., [5]) have assumed a zero thickness “pencil beam” for a DSN antenna. However, on time 
scales less than &10 s, the nonzero diameter of a DSN antenna will modify this pencil beam spectrum. 
The method for calculating this “aperture averaging” effect is presented in Section II, and results are 
given in Section III. A short summary is given in Section IV. 


II. Calculation of Aperture Averaging 

Both theoretical arguments [6] and observational data [2,4,7] support the idea that tropospheric re- 
fractivity fluctuations obey a Kolmogorov power law. The refractivity structure function D n (r) on scales 
up to at least a few tens of kilometers is 


D„{r ) = ([N( x + r) - N(x )] 2 ) = C 2 r 2/3 (1) 

where N(x) is the refractivity at location x (N = n — 1, where n is the index of refraction), and C n 
is the refractivity structure constant. Over time scales up to thousands of seconds, time variations in 
line-of-sight tropospheric delay can be successfully represented by a “frozen flow” model, in which spatial 
variations are convected past the observer by the wind. 

The height dependence of C n is poorly constrained by data. It was modeled by [4] as constant up to 
a 1-km height and zero above that. This “slab height” was subsequently revised to 2 km. 2 Recent WVR 
data from Goldstone [3] show that the observed fluctuations agree with the predictions of this 2-km slab 
model on time scales >400 s but are lower than model predictions on shorter time scales. This result is 
most easily explained by a thickness of the turbulent layer that is significantly larger than 2 km, although 
the data are not adequate to solve for a specific value for this thickness. The effect of a finite thickness 
h of the medium is that fluctuations in the vertical dimension saturate on time scales > h/v w , where v w 
is the wind speed. For time scales At « h/v w , the delay structure function D r oc (At) 5 / 3 , and for time 
scales At » h/v w , D r oc (At) 2 / 3 [4]. The time scales covered in this article are all << h/v w , WVR 
measurements of fluctuations on time scales of 200 s [3] were used to derive the mean C n , once h was 
chosen. For h = 4 km, C n = 3.0 x 10 -8 m” 1 / 3 ; these values were used for all calculations presented in 
this article. Over the course of a year, the structure constant at Goldstone exhibits variations of at least 
a factor of 2 about its mean value [3]. 

The tropospheric delay for a pencil beam looking from location x at elevation angle 9 and azimuth 
AZ is 


1 r°° 

T P encii(x , 0, AZ) = — — / N [x + f(9, AZ, z)\ dz (2) 

sin v J 0 

where r(9,AZ,z) is the vector from the surface to height z in direction (9, AZ). For a DSN antenna of 
diameter d , the delay is averaged over the circular aperture, to give 

4 /*2 7T nd/2 

TDSN (x,6,AZ) = J J £ Tpencii [x + f (£, <f>) , 6, AZ] d£d<j) (3) 

where r(£, (j)) is the vector in the plane of the aperture, starting from the center, of length £ and at angle 
<; f > (the direction chosen for (j) — 0 does not affect the result). Equation (3) assumes that the near- field 


2 R. N. Treuhaft, personal communication, Jet Propulsion Laboratory, Pasadena, California, 1995. 
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DSN beam profile is uniform across its circular cross section. The true profile will be tapered towards 
the edges. As presented in Section III, a change of a factor of two in antenna diameter causes a change 
in the resulting fluctuation level by a factor of «2. Therefore, neglect of beam tapering should cause an 
error of <20 percent in the results reported here. 


A convenient way of calculating and expressing the tropospheric fluctuation spectrum is Allan variance 
(At), defined for a delay process r(t) as [8] 


^(At) = 


(\r(t + 2A t) — 2 r(f + At) + r(£)] 2 ^ 
2 (At) 2 


( 4 ) 


Expanding Eq. (4) and assuming that averaged quantities are independent of time and position gives, for 
the tropospheric delay fluctuations measured by a DSN antenna, 


2 / a _ 3 ( r DSJvW) 4 {roSN(t + At)TDSN(t)) (TDSN(t + 2At)TDSN(t)) 
' ~ ( A<) 2 (At) 2 + : (At) 2 


We can evaluate this expression with Eq. (3) and the frozen-flow assumption. Taking the middle term in 
Eq. (5) as an example, 


16 /*2tt /’27T rd/ 2 rd/ 2 

{TDSN(t "f* A £)tc> 5 jv(£)) — ^2^4 J J J y pencil [ x “b V w At + r(£, </>)] T penc n [x 4" r(£ , (j) )]) 


x d^'d^dfi 


It has been assumed that averaged quantities are position independent, so that 


{^"pencil (^1 » ^5 ■AZ)T penc ii (#2, 0, AZ)) — (Ty enc n(x,6 , AZ)^ ^ ^ [^pencil (*^1 ) AZ) Tp enc n(x2, 0, AZ)} ^ 

= (rfencil&0> AZ ))-\ D T(\xi-X 2 \) 


where D r is the pencil-beam delay structure function, defined as 


D t (t) — ^ I'l'pencili^' A T , 0, AZ) T p ^ nc il(x , 6 , AZ ) ] ^ 


In evaluating Eq. (5), the (r^ encil (x)^ terms add to zero, and ay(At) can be expressed as a four- 
dimensional integral of three D r terms: 
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•ig p2tt p2TT /*cf/2 fd/2 

^ )= «fi -n n « l ^ 


■\d t {A) + 2D t {B)-\d t {C) 


A = \r{t,4>) ~ ?{£! ,<t>')\ 


B = \v w At + r(t,<j>)-TX£',<t>')\ 


> ( 6 ) 


C = |2vtoAt 4- r(£, </>) - f(£', <£')| 


The horizontal distance scales are d and v w A t. As a result, cr y (At) will be suppressed, relative to its 
d = 0 value, for values of At smaller than a few times the wind speed crossing time of t cross = d/v w . 
For a typical wind speed in the lower few kilometers of 8 m/s and for d = 34 m, t cross « 5 s. For 
A t « t cr oss , cry (At) will be much lower than its d = 0 value, due to the averaging of many small-scale 
delay fluctuations over the much larger DSN aperture. 

For any DSN tracking scenario at an observing frequency high enough (>8 GHz) for tropospheric 
delay fluctuations to dominate over those from plasma, the Earth’s troposphere will lie entirely within 
the near field of the antenna. The effective tropospheric volume sampled by the DSN antenna will then 
be cylindrical. For tracking of planetary spacecraft, the antenna will move at a nearly sidereal rate. This 
tracking will cause an effective velocity u e // through the troposphere, at height z and elevation angle #, 
of 


v eff 


Q,z 

sin# 


The Earth’s rotation rate is f2. For z = 2 km (half the turbulent slab height), v e ff = 0.15/ sin# m/s. 
This is much less than the wind speed of 8 m/s used in the calculations. Sidereal tracking, therefore, will 
not significantly modify the results presented in this article. 

The Fresnel length scale l for observations at a wavelength A with phase perturbations occurring at a 
distance L is l « a/A L. For A = 1 cm and L — 2 km (the middle of the turbulent troposphere), l & 4 m. 
On smaller scales, geometric optics is not completely valid. The results for time scales <4 m/s/v w = 0.5 s 
should, therefore, be regarded as approximate. 


III. Results 

Equation (6) was evaluated numerically, for a range of time intervals, and for both 34-m- and 70-m- 
diameter antennas. A wind speed of 8 m/s was assumed. The results for antennas pointed in the zenith 
direction are shown in Fig. 1. Values of cr y (At) for a pencil beam antenna are shown for comparison. 

For time intervals shorter than 10 s, the Allan deviation for measurements made with a DSN antenna 
are less than those made with a pencil beam. For time intervals less than «2 s, a y (At) oc (A £)*% where 
k « 4-1. This proportionality can be roughly understood as follows: The effect of fluctuations over a time 
interval A t are confined to a region of size &v w A t. For v w A t « d, there are N regions ~ ( d/(v w At )) 2 
such regions over the antenna aperture. Averaging the delay fluctuations of these regions will cause a 
reduction by a factor of ~y/N region3 ~ d/(v w At). The slope of the pencil beam a y (At) is « — 1/6, so the 
slope of the DSN cr y (At) will be «+5/6. 

As At decreases, the amount of computer time needed to calculate a y (At) increases rapidly. For small 
At, the terms A, B , and C in Eq. (6) are nearly equal, and high numerical precision is needed in order 
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Fig. 1. The Allan deviation calculated for a pencil 
beam antenna and for antenna diameters of 34 and 
70 m is shown as a function of time interval for 
typical Goldstone conditions and the zenith 
direction. The Allan deviation at other elevation 
angles will be somewhat larger, increasing by a 
factor of 1.8 at a 10-deg elevation angle. The thermal 
noise from a WVR with the characteristics expected 
for the Cassini radio science calibration system is 
included. 


to determine the differences in the D r terms. Furthermore, all derivatives of D n and all but the first 
derivative of D r become singular at zero separation, so that numerical integration methods have difficulty. 

Calculations for At as small as 0.02 s were possible in the zenith direction, because D r (r) could be 
evaluated separately, with an empirically fitted analytical formula used for D r in Eq. (6) to reduce the 
number of dimensions in the integration from 6 to 4. Calculations at other elevation angles were limited 
to a narrower range of time intervals (1-100 s). These calculations showed that a y (At) increases slowly 
with decreasing elevation angle. Relative to its value at the zenith, cr y (At) was 1.16 times larger at a 
30-deg elevation angle and 1.82 times larger at a 10-deg elevation angle. 

All cr y (At) values in Fig. 1 are linearly proportional to the structure constant, C n . Therefore, the 
actual a y (At) values at a DSN site will range up and down from those in Fig. 1 by a factor of at least 2, 
depending upon the season, time of day, and weather conditions. Winter nights tend to have the lowest 
C n values, with summer days having the largest values. 

Tropospheric delay fluctuations at microwave frequencies can be calibrated by WVRs, which measure 
thermal emission from water vapor in the vicinity of its 22-GHz spectral line [9]. However, on short 
time scales, the thermal noise from a WVR becomes the limiting error source. The main vapor-sensing 
channel of current WVRs is at a frequency near 21.5 or 23.8 GHz, where a 1-K brightness temperature 
corresponds to mm jc = 2 x 10“ n s of path delay (c is the speed of light). The thermal path delay 
equivalent noise in a WVR measurement of integration time ti nt for a total system temperature T sys and 
bandwidth BW is 


N(t int ) = 


T S ys 

y/BW ti nt 


2 X 10" 11 S 
IK 


2 x 10- 13 T 100 

s 


( 7 ) 


where Ti 0 o = T sys /100 K, BW 100 = BW / 100 MHz, and t inta = t int /l s. 
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To set a firm lower limit on the effect of WVR thermal noise, we assume that At = ti nt (i.e., we give 
up all information on time scales shorter than At). Future uncooled WVRs currently under design are 
expected to have T sys « 300 K and BW « 400 MHz. 3 The Allan deviation from thermal WVR noise for 
A t 3 = At/1 s (note that <r y (At) = N(t int )y/3/At for thermal noise) then is 


a y (At) « 


5 x 1CT 13 

W 5 


(WVR thermal noise) 


( 8 ) 


This WVR thermal noise curve is plotted in Fig. 1. If we require cr y (At) from WVR thermal noise to be 
at least 3 times smaller than the tropospheric cr y (At) in order to achieve useful calibration, then, for any 
elevation angle >10 deg, the minimum time interval A t m * n for WVR calibration is 


At 


min 


10 s 


Because the slope of the tropospheric cr y (At) is small at At ~ 10 s and the slope of the WVR thermal 
noise <r y (A t) is steep, the minimum useful calibration interval is relatively insensitive to the values of 
T sys and BW. On time scales less than ^10 s, tropospheric fluctuations cannot be calibrated with any 
known technique. 


IV- Summary 

The effect of aperture averaging upon radio science measurements is that tropospheric fluctuations will 
not be important on short time scales. For occultation measurements of Saturn’s rings with Cassini, the 
shortest time scale of interest is the desired physical resolution (wlOO m) divided by the spacecraft orbital 
velocity («10 km/s), or ~0.01 s. 4 These measurements will use an onboard oscillator with cr y (A£) > 
10“ 13 . Even if future flight oscillators were a factor of 10 more stable than the one on Cassini, the 
troposphere would be more stable than the oscillator for At < 1 s. The tropospheric fluctuations on 
At < 10 s appear to form a fundamental lower limit for radio science measurements using an Earth 
antenna, because WVR thermal noise precludes calibration of the fluctuations. 
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